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ABSTRACT 


The  rate  of  production  of  nitric  oxide  in  shock-heated  air  was 
measured  over  the  temperature  range  from  2300°K  to  4500°K.  The  NO 
concentration  was  determined  from  the  NO  infrared  emission  in  the  5.  3 
micron  fundamental  vibration-rotation  band.  The  main  production  of  nitric 
oxide  was  found  to  be  by  the  dissociation  reaction  02  +  M-^420  +  M,  followed 
by  the  atomic  shuttle  reactions  O  +  nL2  NO  +  N  and  N  +  Oz  —  NO  +  O.  The 
overall  rate  constant  k  =  k2«k^  for  NO  production  was  found  to  be  0.40  ± 

0.  12  of  the  rate  constant  compiled  by  Wray. 

The  rate  constant  for  the  bimolecular  reaction  N2  +  02  —  2NO  was 
found  to  be  less  than  one-tenth  the  value  suggested  by  Freedman  and  Daiber. 
In  order  to  explain  the  disagreement  in  the  rate  constant,  experiments  were 
performed  in  shock-heated  nitric  oxide  where  tne  depletion  of  nitric  oxide 
and  the  reaction  products  were  observed.  It  was  found  that  the  NO  depletion 
occurs  mainly  via  the  reaction  2NO  —  NzO  +  O. 

An  analysis  of  the  vibrational  relaxation  of  N0-02-N2  mixtures  is 
presented  and  it  is  shown  that  NO  vibration  lag  is  small  and  does  not  affect 
the  data  analysis.  We  determined  that  the  integrated  band  intensity  of  the 
NO  fundamental  is  124  +  22  atm'  1  cm'2  at  S.  T.  P. 

It  is  also  shown  that  corrections  due  to  the  shock  tube  boundary  layer 
must  be  included  in  order  to  obtain  detailed  agreement  between  experiment 
and  calculations.  The  procedure  for  correcting  chemical  measurements  for 
shock  tube  boundary  layer  growth  is  outlined. 
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I.  INTRODUCTION 


Interest  in  the  chemical  reaction  rates  in  high  temperature  air 

arises  in  many  combustion  studies,  such  as  in  controlled  flame  reactions 

and  in  rocket  motors.  Also,  chemical  changes  in  shocked  air  become 

significant  in  high  speed  flight  above  Mach  6.  Nitric  oxide  (in  equilibrium) 

is  an  important  constituent  of  air  above  2000°K.  ^  Several  authors 

have  tabulated  the  various  reaction  paths  for  the  production  of  NO  starting 

with  N2-02  mixtures.  For  many  of  these  reactions,  the  rate  constants 

were  measured  over  a  small  temperature  interval,  and  theoretically 

extrapolated  over  a  wide  temperature  range.  The  purpose  of  this  work  was 

to  determine  the  rate  of  production  of  NO  in  shock  heated  air.  This  is  an 

integral  type  experiment  in  that  many  reactions  take  place  simultaneously. 

Our  results  were  then  compared  to  calculated  NO  production  rates  in  air 

4 

based  on  the  reaction  rate  constants  compiled  by  Wray. 

An  extensive  high  temperature  shock  tube  study  of  NO  formation  in 

N2-Oz  mixtures,  diluted  in  Ar,  was  performed  by  Glick,  Klein  and  Squire5 

employing  the  single  pulse  shock  tube  method.  The  gas  resided  at  a  high 

temperature  for  several  milliseconds,  and  then  cooled  quickly.  Then  a  gas 

sample  was  taken  from  the  shock  tube  and  analyzed  for  NO.  The  argon  acted 

as  an  additional  catalyst  for  quickly  bringing  the  atomic  oxygen  concentration 

to  equilibrium.  6  Thus,  they  were  able  to  determine  the  rate  constants  for 

the  NO  producing  reaction,  O  +  N2  NO  +  N.  This  reaction,  followed  by 

the  faster  reaction  N  +  02  - *  NO  +  O,  forms  a  pair  of  NO  producing  reactions 

with  the  atoms  acting  as  catalysts.  In  the  case  of  air,  the  atomic  oxygen  is 
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produced  either  by  02  +  M  -*  20  +  M  (where  M  ~  N2  or  02)  or  N2  +  02 
N20  +  O.  The  rate  constant  for  the  latter  reaction  has  not  been  determined 
And  its  importance  for  temperatures  above  2000°K  is  not  known.  Before 
these  atom  reactions  become  important,  the  initial  production  oT  NO  is 
reported  to  be  by  the  reaction  N2  +  02  — *  2NO;  however,  in  the  present 

experiment  this  reaction  is  negligible.  Our  measurements  of  the  rate  of 


fl  >  ^ 

m 


NO  formation  in  shock  heated  air  were  compared  to  calculations  based  on  the 

rate  constants  for  the  above  reactions. 

In  addition  to  the  single  pulsed  shock  tube  method  discussed  above, 
several  other  methods  for  determining  NO  in  shock  heated  air  exist.  The 
ultraviolet  absorption9,  lland  the  NO  infrared  emission  in  the  vibration- 
rotation  band12,  ^techniques  have  been  employed  for  observing  the  NO 
vibrational  relaxation  and  decomposition  rates.  At  lower  temperatures, 

NO  reactions  have  been  studied  in  flames ^and  flow  tubes.  In  this  paper, 
the  NO  infrared  emission  technique  was  employed. 

The  theory  of  the  experiment  is  discussed  in  the  next  section.  The 
experimental  measurements  of  the  NO  infrared  emission  from  shock  heated 
air  follows.  In  the  final  section,  our  measurements  are  compared  to  NO 
production  calculations.  The  important  conclusions  are  summarized. 

There  are  three  appendices  in  this  report.  In  the  first,  the  theory 
oi  vibrational  relaxation  in  a  mixture  of  diatomic  molecules  is  reviewed. 
Numerical  results  are  presented  for  the  vibrational  relaxation  of  a  trace 
amount  of  NO  in  air  over  the  temperature  range  from  2000°K  to  6500  K. 

In  the  second  appendix,  band  intensity  formulation  is  developed  and  applied 
to  the  calibration  of  the  infrared  system.  These  formulas  were  used  in 
Section  III  to  determine  the  NO  fundamental  band  strength.  It  is  also  shown 
that  the  weak  line  approximation  is  valid  for  our  shock  tube  data.  In  the 
final  appendix,  the  procedure  for  correcting  the  chemical  measurements  for 
the  effects  of  shock  tube  boundary  layer  growth  is  developed  and  then  applied 
to  the  NO  formation  chemistry. 
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II.  THEORY  OF  EXPERIMENT 


The  main  purpose  of  this  paper  is  to  determine  the  rate  of  formation 
of  NO  in  shock  heated  air.  The  NO  concentration  was  determined  from  the 
NO  infrared  emission  in  thi  5.  3  micrcn  fundamental  vibration  rotation  band. 

In  order  to  evaluate  the  experimental  data,  the  gas  conditions  downstream  of 
the  shock  and  the  sensitivity  of  the  infrared  system  to  the  NO  concentration 
and  vibrational  temperature  were  determined. 

Experimental  data  have  been  obtained  over  the  shock  velocity  range 
from  2.  0  to  4.  5  mm  per  microsecond,  and  at  initial  air  pressures  from  1  to 
50  torr.  The  shock  heated  gas  temperature  and  density  are  shown  in  Figs.l 
and  2  for  three  chemical  states  of  air  —  frozen  vibration,  equilibrium  vi¬ 
bration  and  frozen  chemistry,  and  full  chemical  equilibrium.  The  equilibrium 
species  concentrations  are  presented  as  a  function  of  incident  shock  velocity 
in  Fig.  3.  These  conditions  were  computed  by  combining  the  Rankine- 
Hugoniot  shock  relations  with  the  equation  of  state  for  air.  The  calculations 
of  thermodynamic  properties  and  composition  of  equilibrium  air  were  based 
on  published  spectroscopic  data.  *5*16  jsj0te  that  in  equilibrium,  oxygen  disso¬ 
ciation  becomes  appreciable  and  the  NO  concentration  reaches  several  percent. 
The  heavy  lines  in  Figs.  1  and  2  indicate  the  range  of  experimental  data. 

The  intensity  of  the  radiation  from  the  entire  5.  3  micron  band  is 
directly  related  to  the  vibrational  energy  even  for  nonequilibrium  vibrational 
distributions.  This  can  be  shown  as  follows:  the  energy  E  invested  in 
vibration  is 

E  =  N  E  (1) 

*  mm 
m 

Nm  and  Em  are  the  number  density  and  energy  of  the  m**1  vibrational  level. 

In  the  harmonic  oscillator  approximation,  this  becomes 
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Fig.  1  Temperature  of  shock  heated  air  as  a  ftinction  of  shock  velocity 
for  several  initial  pressures.  The  heavy  lines  indicate  the  range 
of  the  experimertal  data. 
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Fig.  2 


Density  of  shock  heated  air  as  a  function  of  shock  velocity  for 
several  initial  pressures.  The  heavy  lines  indicate  the  range  of 
the  experimental  data. 
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where  E0  is  the  energy  gap  between  equally  spaced  levels.  In  the  same 
approximation,  the  infrared  emission,  I,  is  due  to  the  radiati*. e  transitions 
between  adjacent  levels. 
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where  hp  =  E  and  t  arc  the  photon  energy  and  the  lifetime  of  the 
.mo  m 

mth  level.  The  lifetime,  r  ,  determined  from  the  dipole  selection  rule 

m 

for  harmonic  oscillators,  is  r  =  T,/m.  (In  this  paper,  we  find  that  t. 

ml  i 

is  approximately  0.  I  second  for  the  NO  fundamental.)  Substituting  for  hp^ 

and  r  in  Eq.  (3),  and  noting  the  similarity  with  Eq.  ,  (2), 
m 


I  = 


/  ,  N  mE 

Am#  m  O 


(4) 


Thus,  in  the  harmonic  oscillator  approximation,  the  infrared  emission  is 
directly  proportional  to  the  vibrational  energy.  In  practice,  the  coupling 
between  vibrational  and  rotational  modes  must  be  included;  however,  the 
emission  from  the  entire  band  remains  within  a  few  percent  of  that  given  by 
Eq.  (4). 

NO  Vibrational  Relaxation  in  Air 

In  the  harmonic  oscillator  approximation,  the  vibrational  energy  is 
related  to  the  vibrational  temperature,  Tv  ,  by  the  l  elation 


ME 

E “  "E  /kT 
o'  v 
e 

where  N  is  the  number  density  of  NO  molecules  and  k  is  the  Boltzmann 
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constant.  It  has  been  shown  that,  in  an  abrupt  change  in  the  gas  tempera¬ 
ture  causing  vibrational  lag,  the  vibrational  level  population  adjusts  so 
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that  a  vibrational  temperature  Tv  exists  during  the  relaxation 
process. 

In  order  to  relate  the  infrared  emission,  I,  to  the  NO  concentration, 

N,  by  Eq.  (5),  the  NO  vibrational  temperature  must  be  known.  In  a  chemical 
experiment,  the  infrared  signal  can  be  used  as  a  direct  measure  of  the  NO 
concentration  only  when  the  NO  vibrational  and  translational  temperatures 
are  equal.  However,  immediately  behind  the  shock,  the  vibrational 
temperature  of  NO  (also  of  N2  and  02)  is  less  than  he  kinetic  temperature, 
and  the  infrared  emission  is  less  than  that  for  equilibrium.  Thus,  a  correct 
interpretation  of  the  infrared  signal  requires  a  knowledge  of  the  vibrational 
relaxation  rates.  The  general  problem  of  vibrational  coupling  between  NO, 
O2  and  N2  is  considered  in  detail  in  Appendix  A.  The  results  of  these  calcu¬ 
lations  are  reviewed  in  Section  IV  where  they  are  applied  to  the  experimental 
data.  Fortunately,  the  vibrational  lag  correction  to  the  data  is  small. 

Chemical  Reactions  for  NO  Formation 

The  chemical  reactions  leading  to  NO  production  have  been  reviewed 
in  the  literature  2-4  while  over  a  score  of  possible  reactions  are  tabulated, 
there  are  only  a  few  important  reactions  and  these  can  be  arranged  in  three 
groups: 

(1)  Direct  production  of  NO  by  the  reaction 

Reaction  1  N2  +  O 2  “*  2NO  (6) 

(2)  Atomic  shuttle  reactions  which  convert  N2  and  O2  to  2NO 

Reaction  2  O  +  N2  -  NO  +  N  (7) 

Reaction  3  N  +  02  -*  NO  +  O  (8) 

These  reactions  have  an  incubation  time  until  there  is  sufficient  production 
of  atomic  oxygen;  at  such  time  they  become  the  dominant  reactions  producing 
NO. 

(3)  Atom  producing  reactions 


Reaction  4 

Oz  +  M  —  20  +  M 

(9) 

Reaction  3 

N2  +  M  —  2N  +  M 

(10) 

Reaction  6 

NO  +  M  —  N+O  +  M 

(ID 

M  can  be  any  of  the  atoms  or  molecules  shown  in  Reactions  1  to  6. 
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The  reaction  N2  +  O2  —  N^O  +  O  could  also  be  a  source  of  atomic  oxygen, 
but  its  rate  is  not  known. 

Our  experimental  measurements  were  compared  to  calculated  NO 
formation  rates  which  were  based  on  the  rate  constants  compiled  by  Wray 
and  shown  in  Table  I.  These  rates  were  varied  in  order  to  determine 
the  NO  production  sensitivity  to  the  rate  constants.  A  comparison  of  these 
calculations  with  the  data  is  presented  in  Section  IV. 
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TABLE  1  RATE  CONSTANTS  FOR  CHEMICAL  PROCESSES 


III.  EXPE1  1ENTAL  MEASUREMENTS 


# 


Experimental  Apparatus 

The  experiments  were  carried  out  in  a  1.  5  inch  diameter  stainless 
steel  shock  tube  composed  of  a  driver,  a  22  foot  long  low  pressure  section, 
and  a  large  expansion  tank  at  the  upstream  end.  X-scribcd  stainless  steel 
diaphragms  were  ruptured  . jy  high  pressure  H2  or  H2-N2  mixtures.  A  10 
mil  cellophane  diaphragm  separated  the  low  pressure  section  from  the  expansion 
tank.  To  maintain  a  low  impurity  level  and  also  to  reduce  the  pumping  time, 
special  valves  were  developed  which  permitted  changing  diaphragms  at  the 
driver  and  expansion  tank  while  the  shock  tube  was  under  a  partial  vacuum. 
Thus,  room  air  did  not  enter  the  shock  tube  between  runs.  These  valves 
had  an  auxiliary  chamber  in  which  new  diaphragms  could  be  pumped  on  and 
out-gassed  before  insertion  in  the  shock  tube.  The  shock  tube  was  evacuated 
to  a  pressure  of  2  x  10“6  torr  or  lower  by  an  oil  diffusion  pump  with  a  liquid 
N2  cooled  trap.  The  leak  rate  after  sealing  off  the  shock  tube  was  approxi¬ 
mately  0.  3  microns  Hg  per  minute.  Since  the  inital  gas  was  dried  room 
air,  the  effects  of  this  small  air  leak  on  the  chemistry  was  negligible. 

The  shock  velocity  was  measured  with  a  series  of  platinum  resistance 
heat  transfer  gauges  spaced  several  feet  apart.  The  heat  transfer  signals 
were  displayed  on  a  raster  oscilloscope.  The  distance  between  gauges  was 
measured  to  t  0.01  inch  and  the  time  between  stations  was  determined  to 
1/2  microsecond.  The  maximum  error  in  the  shock  velocity  determination 
was  +  0.65%.  The  shock  velocity  attenuation  was  constant  along  the  length 
of  the  tube  and  was  only  a  function  of  shock  velocity:  0.  20  mm  per  micro¬ 
second  per  meter  length  of  tube  at  the  high  velocities,  and  0.  074  mm  per 
microsecond  per  meter  length  of  the  tube  at  the  low  velocities.  Where  the 
effects  of  shock  attenuation  became  important,  the  data  were  not  analyzed  in 
detail. 

A  schematic  diagram  of  the  section  of  the  shock  tube  for  the  infrared 
measurements  is  shown  in  Fig.  4.  This  section,  located  18  feet  from  the 
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experimental  arrangement 


\ 

\ 


I 
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4  Schematic  diagram  of  shock  tube  section  with  infrared  detection 
systems. 


diaphragm,  was  made  of  lucite,  a  poor  infrared  reflector.  At  b  microns 
wavelength,  lucite  is  opaque  with  less  than  5%  reflectivity.  In  the  initial 
runs,  metal  walls  were  used,  but  they  gave  extraneous  signals  due  to  reflec¬ 
ted  radiation  originating  downstream  of  the  window.  Two  infrared  optical 
systems  we.-e  located  in  the  same  plane  normal  to  the  shock  tube  axis.  The 
detector  in  each  system  was  imaged  on  the  inner  surface  of  the  CaF£  infra¬ 
red  window  with  an  8  inch  diameter,  10  inch  focal  length  spherical  mirror 
set  at  approximately  unit  magnification.  An  infrared  filter  was  placed  in 
front  of  each  detector  to  isolate  the  band.  The  CaF2  infrared  windows  were 
radiused  to  the  inside  diameter  of  the  shock  tube.  System  I,  used  for  good 
spatial  resolution,  had  a  1  inch  by  8  inch  lucite  aperture  over  the  mirror, 
the  narrow  dimension  being  along  the  axis  of  the  shock  tube.  The  1  mm 
by  5  mm  active  area  of  the  (Philco)  InSb  infrared  detector  defined  the 
image  size  with  the  1  mm  width  imaged  along  the  axis  of  the  shock  wave. 
There  were  no  slits  at  the  shock  tube.  Measurements  with  this  InSb  detec¬ 
tor  required  rise  time  corrections  when  the  signal  rise  time  was  faster 
than  four  microseconds.  For  the  low  shock  velocity  runs,  System  II  used 
the  entire  mirror  area  and  a  5  mm  diameter  (Philco)  InSb'Infrared  detector. 
This  system  has  the  maximum  sensitivity  and  was  used  to  investigate  the 
impurity  level  and  the  initial  NO  low  level  signals.  For  the  high  shock 
velocities,  the  aperture  on  the  mirror  was  reduced  to  1"  by  8"  and  a 
(Westinghouse)  Ge:Au  detector  with  a  1mm  by  3mm  effective  sensitive 
area  and  with  0.  7  microsecond  overall  system  rrse  time  was  used.  The 
outputs  of  the  infrared  detectors  were  fed  directly  to  standard  low  noise 
high  gain  electronics  and  then  displayed  on  oscilloscopes. 

Figure  5  shows  the  spectral  response  of  the  InSb  and  Ge:Au  detectors 
and  the  interference  filters  used  in  the  experiment. 

In  most  of  the  shock  tube  runs,  the  initial  gas  was  room  air  passed 
through  an  ascarite  filter  and  a  liquid  nitrogen  trap  to  remove  I^O  and 
CC^.  A  series  of  check  runs  in  synthetic  air  made  from  Matheson  high 
purity  N2  and  O2  ehowed  no  difference  in  results  when  dried  room  air 
was  used.  The  H2  used  in  the  driver  was  replaced  by  He  and  a  series 
of  runs  were  repeated.  The  purpose  of  these  runs  was  to  ensure  that 
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RELATIVE  RESPONSE  (PERCENT) 


Fig.  5  The  relative  response  for  indium  antimonide  and  gold-doped 

germanium  infrared  detectors  as  a  function  of  wavelength.  The 
relative  transmittance  of  the  filters  used  in  the  experiment  is 
also  shown  as  a  function  of  wavelength.  The  curve  for  the  gold- 
doped  germanium  detector  was  supplied  by  the  manufacturer;  the 
other  curves  were  obtained  from  laboratory  measurements. 


hydrogen  compounds  from  the  previous  run,  that  might  be  produced  in 
the  mixing  region  and  deposited  on  the  walls  of  the  shock  tube,  would  not 
contribute  to  the  infrared  emission.  There  were  no  detectable  differences 
using  the  H 2  or  He  driver.  Another  series  of  runs  was  also  made  in  pure 
N2  and  O2  to  check  for  infrared  emission  from  impurities  in  these  gases, 
and  from  boundary  layer  impurities  6uch  as  lucite  ablation  products.  Fcr 
the  and  runs,  only  System  II  with  the  5  mm  diameter  InSb  detector 
showed  a  small  signal,  while  System  I,  eight  times  less  sensitive  had  no 
indication  above  noise.  The  signals  remained  relatively  constant  with 
shock  velocity  using  either  N2  or  O^.  This  radiation  might  have  been 
due  to  boundary  layer  effects  or  impurities  in  the  gases;  the  intensity 
was  the  equivalent  to  NO  radiation  from  a  2  to  4  x  10  *  NO  concentra¬ 
tion.  The6e  signals  could  be  increased  by  increasing  the  shock  tube  base 
pressure  to  several  microns  Hg  before  filling  the  shock  tube.  This  in¬ 
creased  signal  was  consistent  with  about  10  ppm  of  water  vapor. 

/ 

Infrared  System  Calibration 

The  output  signal  from  the  infrared  system  was  calibrated  to  give 
the  NO  concentration  in  the  shock  tube.  In  a  series  of  shock  tube  runs  using 
NO-Ar  mixtures  the  NO  infrared  signal  was  determined  as  a  function  of  the 
NO  temperature,  concentration  and  partial  pressure.  Figure  6  shows  data 
oscillograms  for  1,  10  and  100  percent  NO  in  Ar  runs.  The  voltage  of  the 
level  portion  (after  either  NO  vibrational  relaxation  or  system  response 
time)  normalized  to  l  torr  NO  partial  density  is  plotted  as  a  function  of 
the  temperature  in  Fig.  7. 

The  curve  in  Fig.  7  is  the  calculated  NO  emission  using  Eq.  (B-12) 
which  assumes  an  optically  thin  gas;  that  is,  the  contribution  from  each 
molecule  is  additive.  The  height  of  the  curve  was  adjusted  for  the  best  fit 
to  the  data.  Note  the  excellent  agreement  with  the  computed  curve  except 
for  the  runs  with  100%  NO  at  an  initial  pressure  of  5  torr,  which  are  about 
30%  low.  This  decrease  is  due  to  self-absorption  effects  in  which  the  center 
of  the  lines  become  black.  For  the  shock  tube  runs  in  air,  the  partial  pres¬ 
sure  of  nitric  oxide  behind  the  shock  is  less  than  1/50  this  value.  In  addition, 
there  is  pressure  broadening  by  O2  and  N£.  Thus,  the  air  cata  does  not 
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NORMALIZED  INFRARED  OUTPUT  -  mv 


Fig.  7  NO  infrared  emission  as  a  function  of  temperature.  The  solid 
curve  shows  the  computed  intensity  variation  with  temperature 
assuming  the  weak  line  approximation.  The  data  was  normalized 
at  (1/760)  STP. 
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require  corrections  for  self-absorption;  the  infrared  emission  is  propor¬ 
tional  to  NO  concentration  with  the  temperature  dependence  given  in  Fig.  7, 
see  also  Appendix  B. 

The  infrared  systems  were  also  calibrated  and  checked  periodically 
with  a  black  body  light  source.  The  primary  purpose  of  the  calibration  was 
to  ensure  that  there  were  no  changes  in  the  sensitivity  during  the  experi¬ 
ment.  A  portion  of  the  shock  tube  was  ’•eplaced  with  a  black  body  at  500°C. 
A  chopper  wheel  in  front  of  the  source  furnished  an  AC  signal.  In  order  to 
vary  the  energy  incident  on  the  infrared  detector,  apertures  were  placed 
over  the  focusing  mirrors,  always  keeping  the  black  body  temperature  at 
500OC.  This  checked  the  linearity  of  the  systems,  and  furnished  signals 
with  the  same  range  of  intensity  as  in  the  shock  tube  runs. 

Figure  8  shows  the  black  body  calibration  signal  as  a  function  of  the 
exposed  fraction,  X,  of  the  mirror  area.  The  data  reproduced  to  within 
10%  over  the  experimental  period.  This  calibration  is  for  System  I  with  the 
xnSb  cell  and  filter  B.  Similar  calibrations  were  obtained  for  the  other 
systems. 

NO  Fundamental  Band  Intensity 

The  band  intensity  of  the  5.  3  micron  NO  fundamental  band  can  be 
determined  from  the  equilibrium  NO  emission  measurements  presented  in 
Fig.  7.  The  previous  discussions  showed  that  the  gas  was  optically  thin. 
Thus,,  ve  can  use  the  formulas  developed  in  Appendix  B.  We  rtart  with 
Eq.  \-3-12)  which  gives  the  band  emission  coefficient  cem  as  a  function  of 
the  NO  emission,  S,  and  the  black  body  emission,  Sg, 


c 

om 


G)®«{ 


CT  T  Ib  dv 


(B- 12), 


(12) 


Numerical  values  of  the  spectral  and  temperature  dependence  for 
the  bracketed  factor  are  presented  in  Fig.  B-l.  Use  of  the  ratio  S/Sg 
eliminates  the  need  to  determine  the  optical  solid  angle  and  the  absolute 
magnitude  of  the  filter  transmission,  r,  and  detector  sensitivity,  o  . 


INFRAREO  SIGNAL 


FRACTION  OF  MIRROR  AREA  EXPOSED,  X 

Fig.  8  The  static  calibration  of  the  infrared  recording  system  is  accom¬ 
plished  by  replacing  the  shock  tube  with  a  500  C  black  body  and 
varying  the  area  of  the  collecting  mirror.  The  observed  signal  is 
plotted  as  a  function  of  the  fraction  of  the  mirror  area  uncovered. 
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To  illustrate  how  Eq.  (12)  is  used,  we  will  determine  cem  with 
the  data  at  3000°K.  From  the  data  shown  in  Fig.  7,  N  =  3.6  x  10  . 

X  =  1/2  *,  t  =  3.81  cm,  and  S  =  100  mV  at  3000°K.  From  Fig.  8, 

O'  o 

(S  /X)  =  1.45  x  104  mV.  The  value  of  the  bracketed  factor  at  3000  K 

is  9.25  x  10”8,  see  Fig.  B-l.  Inserting  these  values  into  Eq.  (12), 

c  becomes 
em 


(10Z)  (2  ?r)(9.  25  x  10  ) 

(3.81)  (1.45  x  104)  (3.6  x  1016) 


=  2.  96  x  10 


•26  ergs/sec 
part.-ster 


c  was  calculated  using  the  data  in  Fig.  7  by  a  least  squ  ires  fit  to 
em 

Eq.  (12). 


c  =  (2.97  +  0.4)  x  10' 
em 


(statistical  error) 


This  uncertainty  is  due  to  the  scatter  of  the  NO  emission  data.  In 
addition,  we  have  considered  the  systematic  errors  listed  in  Table  II.  The 
overall  uncertainty  in  cem  due  to  systematic  errors  is  15%.  Combining 
the  systematic  errors  and  the  statistical  variance,  the  overall  error 
becomes 


c  =  (2.97  +  0.  54)  x  10 
em 


-26  ergs/sec 

particle-ster 


The  NO  band  intensity  can  be  obtained  directly  from  the  emission 
coefficient.  The  band  intensity  is  defined  as  the  integral  band  absorption  of 
an  optically  thin  gas  normalized  to  1  cm3  of  NO  at  standard  temperature  and 
pressure.  In  terms  of  cem»  the  band  intensity  becomes 


/  Hdi 
band 


N  v  i  c 
o  1  em 


v  =  1876  cm"*,  is  the  wave  number  at  the  center  of  the  band;  h  =  Planck's 
constant,  c  =  velocity  of  light  and  Nq  =  2.69  x  1019  (Loschmidt's  Number). 
Inserting  our  value  for  cem,  the  band  intensity  becomes 
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TABLE  II 


SYSTEMATIC  ERRORS 


Source 

Type 

Value 

c  Uncertainty 

em  7 

Black  Body 

Temperature 

-2.  5°C  at  5C0°C 

J  2% 

Black  Body 

Emissivity 

-  1%  at  99% 

Calibration  Mirror 

Reflectivity 

3% 

<1% 

Filter 

Spectral  Response 

<0.  2% 

<1% 

Filter 

Transmittance 

3% 

4% 

Detector 

Spectral  Response 

<1% 

<12% 

Detector 

Relative  Amplitude 

<5% 

<7% 
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band 


-1  -2 

jjjiv  =  124  ±  22  atm  cm  . 


The  values  of  the  5.  3pt  fundamental  NO  band  intensity  obtained  by  others 
are  shown  in  Table  III.  These  values  fail  into  two  groups  --  (1)  high 
values  and  (2)  low  values.  The  high  values  are  in  general  agreement  with 


the  original  work  of  Havens  who  used  the  optical  dispersion  method,  a 


simple  and  reasonably  accurate  way  of  determining  the  fundamental  band 
intensity  of  a  diatomic  molecule.  The  work  of  Penner  and  Weber  was 
repeated  by  Varanasi  and  Penner^  where  the  latter  obtained  a  higher 


value.  No  explanation  for  the  difference  in  the  two  reported  values  was 

21 

presented.  The  low  value  obtained  by  Vincent-Gcisse  is  difficult  to 
r  -,c  27 

explain.  Fukudaf  and  Breeze  and  Fcrriso  worked  at  intensity  levels 


where  the  line  centers  were  black  and  the  weak  line  approximation  was 
no  longer  valid.  This  would  be  an  important  correction  and  should  have 
been  considered  in  the  reduction  of  their  data. 


NO  Production  in  Shock  Heated  Air 

We  will  now  present  the  measurements  of  NO  production  in  shock 
heated  air.  Only  the  "raw"  uncorrccted  data  is  discussed  in  this 
section.  Corrections  and  analysis  of  the  data  are  presented  in  Section  IV. 

Measurements  of  NO  production  in  shock  heated  air  are  shown  in 
Fig.  9.  The  NO  infrared  emission  is  given  as  a  function  of  laboratory  time 
for  several  shock  velocities  in  the  range  of  2,28  to  4.02  mm/  pi  sec  and  at 
pressures  of  1,  5  and  50  torr.  The  shock  heated  gas  conditions  for  the 
range  of  air  data  are  presented  in  Figs.  1,  2  and  3.  The  i'lfrarcd  signal 
shows  an  initial  incubation  period  with  a  low  voltage  level;  it  then  rises 
rapidly  to  a  peak  value.  At  shock  velocities  below  2.  6  mm/  pi,  sec  the 
infrared  signal  rises  monotonically,  while  at  higher  velocities  the  NO 
concentration  peaks  above  equilibrium  and  approaches  equilibrium  from 
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TABLE  III 


INTENSITY  MEASUREMENTS  FOR  THE  FUNDAMENTAL  BAND  OF  NO 
Values  shown  are  referred  to  density  at  STP 


Author 

Method 

Intensity  cm’^atm  * 

Havens  (19 38) ^ ^ 

0,'tical  Dispersion 

121 

19 

Dinsmore  (1949) 

Zero  Path  Extrapolation 

145  ±  29 

Penner  and 

Weber  (1953) 

Piessure  Broadening 

70  ±  7 

Vincent-Geisse 

(1954)^ 

Line  Shape 

82 

Schurin  and  77 

Clough  (1963)  ^ 

.^elf  Broadening 

111  ±  7 

James  (19  64)  ^  ^ 

Line  Shape 

138  ±  6 

Breeze  and 

Ferriso  (1964) 

Integrated  Emission 

76  ±  7 

Fukuda  (1965)^ 

Integrated  Emission 

70  to  77 

Ford  and 

Shaw  (7.965)  ° 

Integrated  Absorption 

115  ±  12 

Abels  and  ,7 

Shaw  (1966) 

Line  Shape 

122  ±  6 

Varanasi  and 

Penner  (1966) 

Pressure  Broadening 

128  ±  10 

This  Experiment 

Integrated  Emission 

124  ±  22 

-• 
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above.  Figure  10  indicates  features  of  the  data  to  be  analyzed:  (1)  the 
incubation  time,  (2)  the  time  to  reach  1/4,  1/2  and  3/4  the  peak  intensity, 
(3)the  NO  concentration  at  peak  intensity.  Only  qualitative  analysis  was 
performed  on  the  data  after  the  peak  intensity  because  of  the  uncertain 
corrections  due  to  the  effects  of  the  side  wall  boundary  layer  and  shock 
attenuation. 

The  incubation  time  may  be  defined  in  t”  ,  ways:  (1)  the  time  at 
which  a  tangent  to  the  fast  rising  portion  intersects  the  zero  level,  and 
(2)  the  time  at  which  the  infrared  signal  becomes  larger  than  the  back¬ 
ground  level.  Essentially  the  same  incubation  time  is  obtained  using 
either  definition.  Figure  11  shows  the  observed  incubation  (laboratory) 
time  as  a  function  of  incident  shock  velocity.  The  data  were  binary  scaled 
to  an  initial  pressure  of  5  torr  by  multiplying  the  observed  time  by 
(Pj/5),  where  Pj  is  the  initial  pressure  in  torr.  The  data  show  that  the 

incubation  time  has  bimolecular  scaling.  The  solid  curves  show  the  N? 

2Q  .  30  .  ^ 

vibrational  relaxation  time  in  pure  ^  and  in  air.  The  possible 
connections  between  the  incubation  time  and  the  N^  vibrational  relaxation 
time  ire  discussed  in  the  next  section. 

The  NO  mole  fraction  at  the  incubation  time  as  determined  from  the 
infrared  signal  is  shown  in  Fig.  12.  As  mentioned  above,  at  this  time  the 

infrared  signal  becomes  larger  than  the  background  signal,  corresponding  to 

-4 

a  NO  mole  fraction  of  2  to  4  x  10  .  Also  shown  in  Fig.  12  is  the  predicted 

mole  fraction  of  NO  using  the  reaction  mechanism  and  rate  constants  given 
in  Table  I.  ^  For  shock  velocities  below  2,  4  mm/ p  sec  the  data  is  consistent 
with  the  predicted  NO  production.  However,  by  3.  0  mm/p  sec  there  is  over 
an  order  of  magnitude  difference.  This  difference  in  the  early  NO  production 
is  the  only  major  discrepancy  between  our  experimental  data  and  predictions; 
and  is  due  to  the  predicted  extra  NO  production  by  the  reaction  N^  +  O^-*  ^O. 

After  the  incubation  period  the  infrared  signal  rises  rapidly  to  the 
peak  value.  The  time  for  the  infrared  signal  to  reach  one-quarter  peak 
height(T  i  pk),  half  peak  height  (Tipk)  and  three-quarters  peak  height  (T3/4pjc) 
are  shown  in  Figs.  13,  14  and  15.  The  half  peak  intensity  point  is  nearly 
half  the  peak  NO  mole  fraction.  The  temperature  is  higher  at  the  half  peak 
point  than  at  the  peak,  thus,  the  NO  emission  per  molecule  is  larger.  This 
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10  Schematic  drawing  of  the  air  data  defining  the  various  experimental 
measurements. 


MOLE  FRACTION  OF  NO  AT  INCUBATION  TIME 


Fig.  12 


Computed  and  measured  NO  mole  fraction  at  the  incubation  time 
shown  in  Fig.  1 1. 
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SHOCK  VELOCITY-mm^r-C 

Fig.  13  Time  for  NO  emission  to  reach  1/4  of  the  peak  value.  Data  uncor 
rected  for  electronic  response  time. 
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increase  is  largely  compensated  by  the  lower  density  at  the  higher 
temperature.  The  measured  time  to  reach  the  peak  intensity  is  not  pre¬ 
sented  because  of  the  large  uncertainty  in  estimating  thir>  time. 

The  peak  infrared  emission,  presented  in  Fig.  16t  can  be  corre¬ 
lated  with  the  NO  concentration.  The  NO  emission  calibration  as  a  function 
of  temperature  is  given  in  Fig.  7.  The  gas  temperature  at  peak  emission 
was  not  measured,  but  was  obtained  from  computer  calculations.  Figure  17 
shows  the  mole  fraction  of  NO  at  the  peak  intensity.  The  curve  shows  the 
calculated  peak  value.  This  curve  is  independent  of  initial  pressure  for 
Pj  S  50  torr,  since  3-body  reactions  make  a  negligible  contribution.  The 
peak  NO  mole  fraction  is  greater  than  that  for  chemical  equilibrium  (see 
Fig.  3),  an  effect  predicted  by  Duff  and  Davidson,  ^  which  is  due  to  local 
equilibrium  between  the  molecules  O^,  N^  and  NO  through  Reactions  ) 
and  (3),  while  O 2  dissociation  proceeds  at  a  slower  rate  through  Reaction  (4). 
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NO  MOLE  FRACTION 


Fig.  i7  The  NO  mole  fraction  at  the  peak  intensity  compared  to  calculation 
using  the  reactions  and  rate  constants  of  Table  I. 
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IV.  DISCUSSION  OF  RESULTS 


The  experimental  determination  of  the  NO  production  rate  in  shock 
heated  air  will  now  be  compared  with  computer  calculations  using  the 
reaction  mechanisms  and  chemical  rates  shown  in  Table  I.  The  major 
difference  between  the  experimental  data  and  computer  calculations  it 
during  the  incubation  period  before  the  atom  shuttle  reactions  become 
important,  the  measured  NO  radiation  being  much  less  than  predicted, 
see  Fig.  12. 

Note  that  the  incubation  time  occurs  at  approximately  the  ^  vibra¬ 
tion  relaxation  time.  While  this  might  be  a  coincidence,  it  does  suggest 
several  reasons  for  the  discrepancy.  We  have  examined  three  possible 
reasons  for  this  lack  of  agreement:  (1)  NO  vibration  lag  reducing  the  infra¬ 
red  emission  at  early  times;  (2)  wrong  reaction  rates  in  Table  I,  in 
particular  Reaction  (1)  N£  ■'  O2  —  2NO;  and  (3)  coupling  of  vibration  and 
chemistry  in  the  reactions  ^  -*  2NO  and  O  +  ^  -*  NO  +  O.  The 

coupling  effect  assumes  that  the  chemistry  cannot  proceed  until  N2 
vibration  approaches  equilibrium  conditions.  First  we  will  consider  these 
three  items  and  then  return  to  a  detailed  comparison  of  data  with 
calculations.  It  will  be  shown  that  NO  vibration  lag  is  unimpor;ant,  and 
that  Reaction  (1)  N2  +  O^  —  2NO  in  Table  I  should  be  deleted.  It  is 
difficult  to  ascertain  whether  vibration  and  chemistry  are  coupled  from 
the  present  data. 

Vibrational  Relaxation  of  NO 

The  NO  concentration  is  determined  from  the  3  micron  band 
intensity  assuming  the  NO  vibrational  and  translational  temperatures  are 
equal.  To  evaluate  the  limitations  of  this  assumption,  we  must  know  the 
NO  vibrational  relaxation  rates.  During  the  incubation  period,  the  NO 
concentration  is  low,  thus  the  NO  vibrational  relaxation  depends  on 
collisions  with  O2  and  N£»  If  the  NO-N2  vibration-vibration  exchange 
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were  last,  NO  would  follow  the  N2  vibrational  temperature.  Since  the 
NO  infrared  emission  per  molecule  is  an  increasing  function  of  the 
vibrational  temperature,  the  radiation  could  be  much  less  during  the 
incuoation  period,  until  the  N2  vibrational  relaxation  is  complete. 

The  general  problem  of  the  vibrational  excitation  of  a  small  amount 
of  NO  in  02-N2  mixtures  is  discussed  in  detail  *n  Appendix  A.  For  trace 
amounts  of  NO,  NO  is  vibrationally  excited  by  thermal  collisions  with  N^ 
via  the  reaction  N2  +  NO  —  N2  +  NO*.31  In  addition,  vibrational  energy 
transfer  to  N2  via  the  reaction  NO*  N2  —  NO  +  is  slow  enough,31  so 
that  the  NO  vibrational  temperature  was  above  that  for  N2  during  the  relaxa¬ 
tion  process.  Figure  11  shows  the  calculated  effects  of  the  NO  vibrational 
relaxation  on  the  NO  infrared  emission.  The  dashed  curves  give  the  times  at 
which  the  ratio  of  the  infrared  emission  to  the  vibration  equilibrium  emission 
is  0.60,  0.70,  0.80  and  0.90.  These  curves  were  determined  using  the 
vibrational  energy  relaxation  curves  in  Figs.  A2-A4.  Also  shown  in  Fig.  11 
are  the  incubation  times.  Note  that  at  the  incubation  time  the  correction  to 
the  infrared  emission  due  to  vibration  lag  is  less  than  10%  over  the  shock 
velocity  range  from  2.0  to  2.8  mm  per  microsecond.  Thus,  NO  vibration 
lag  does  not  account  for  the  low  infrared  emission  during  the  incubation 
ptriod. 

Wrong  Rate  Constant  for  N2  +  02  —  2NO 

A  second  possible  reason  for  the  lack  of  agreement  during  the 
incubation  period  is  that  the  rate  constant  used  for  the  reaction  N?  +  0?  — 

2NO  is  incorrect.  This  reaction  produces  most  of  the  NO  during  the  incuba¬ 
tion  time.  If  the  rate  were  at  least  an  order  of  magnitude  smaller,  then 
there  would  not  be  any  disagreement.  The  rate  constant  used  was  suggested 
by  Freedman  and  Daiber,  11  and  Wray  and  Teare^  who  studied  the  reverse 
reaction  and  measured  the  decrease  of  NO  in  shock  heated  NO-Ar  mixtures. 
They  proposed  the  reaction  2NO  —  N2  +  02  for  the  NO  depletion;  however, 
the  reaction  products  were  not_  observed.  The  alternate  reaction1*1 
-NO  —  N20  +  O  could  also  account  for  the  depletion  of  NO  in  shock  healed 
NO.  If  this  latter  reaction  were  dominant,  then  the  rate  constant  for 
Reaction  (1)  would  be  uncertain. 
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Several  measurements  on  shock  heated  NO  were  made  in  order  to 

check  the  reaction  mechanisms  proposed  by  Freedman  and  Daiber**  and 

9 

Wray  and  T  "iare.  In  addition  to  measurements  of  the  NO  concentration, 

we  observed  atomic  O  and  N^O.  The  NO  and  N^O  concentrations 

were  determined  by  the  emission  in  the  5.  3  and  4.  6  micron  bands, 

respectively.  Atomic  O  was  determined  from  the  chemiluminescent 

reaction  NO  +  C  *  N02  +  (hv  ),  by  observing  radiation  in  a  band  from  0.  5 

to  0.  6  microns.  The  initial  O  production  is  directly  proportional  to  this 

chemiluminescent  radiation,  since  the  initial  NO  concentration  does  not 

vary.  Figure  18  shows  some  of  the  measurements  of  shock  heated  NO 

and  NO-Ar  mixtures.  Note  that  N,0  and  O  increase  linearly  with  time, 

indicating  directly  che  reaction  NO  +  NO  —  N20  +  O.  Figure  19  shows  the 

results  of  analyzing  the  initial  decay  of  NO  and  the  initial  rise  of  O.  The 

large  error  bars  on  the  O  data  are  due  to  incomplete  calibration. 

Note  that  both  rates  are  essentially  the  same.  The  rate  constant 

for  the  depletion  of  NO  is  equal  to  that  for  the  increase  in  O,  and  also 

equal  to  that  suggested  by  Freedman  and  Daiber,  11  the  solid  curve. 

No  observation  of  the  products  of  the  reaction  2NO  —  N2  +  02  was 

made.  The  rate  constant  for  this  reaction  is  unknown,  because  it  was  not 

11  Q 

measured  by  Freedman  and  Daiber,  and  Wray  and  Teare.  The  initial 
slope  of  our  shocked  air  data  can  be  used  to  obtain  an  upper  limit  to  this 
reaction  which  at  4000°K  is  <  0.  10  the  rate  given  in  Table  I.  For  the 
purpose  of  air  chemistry  this  rate  constant  is  small  enough  so  that 
Reaction  (1)  can  be  deleted. 

Coupling  of  Vibration  and  Dissociation 

A  third  reason  for  the  discrepancy  between  prediction  and  observed 
initial  NO  production  may  be  due  to  coupling  between  vibration  and  chemistry. 
The  rate  constants  for  Reactions  (1)  and  (2)  could  depend  on  the  vibrational 
temperature,  and  the  expected  rate  constant  would  not  be  attained  until  N? 
vibrational  equilibrium.  Thu6,  Table  I  gives  the  rate  constants  for  the 
reactions  02~  +  N^  —  2NO  and  O  +  N^  -  NO  +  N,  where  N^  'indicates 
equilibrium  vibration.  Previously  we  showed  that  the  reaction  O,  +  N,  — 

2NO  can  be  ignored;  thus,  we  must  investigate  only  the  latter  reaction  for 
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Fig.  18  Measurements  of  O  and  N2O  production  in  shock  heated  NO. 

O  determined  by  O  +  NO  —  NO2  +  h  v  chemiluminescent  reaction, 
NO  and  N2O  concentrations  by  the  5,3  and  4.6  micron  emission, 
respectively. 
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Fig.  19  Experimentally  determined  rate  constant  for  the  initial  dissociation 
of  NO,  open  circles,  and  production  of  atomic  O,  vertical  bars. 
Curve  gives  Freedman  and  Daiber  measured  rate  constant  for  NO 
depletion. 
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the  effects  of  N,  vibration  lag.  The  coupling  of  vibi  a*  on  and  chemistry 
has  been  observed  in  dissociation  and  a  similar  connection  might  occur 

here.  The  best  place  to  check  for  the  effect  of  the  coupling  is  at  the  end  of 
the  incubation  period  l>y  comparing  the  NO  intensity  at  that  time  to  the  pre¬ 
dicted  value.  Since  the  observed  and  predicted  NO  intensity  is  small  at  the 
incubation  time,  we  cannot  obtain  information  on  vibration-dissociation 
coupling  with  the  present  data. 

It  is  interesting  to  note  that  N^  has  the  slowest  molecular  vibration 
relaxation  rate;  thus,  this  reaction  is  ideally  suited  for  studying  vibration 
coupling  to  chemistry.  The  possibility  of  such  coupling  in  Reaction  (2) 
could  be  investigated  experimentally,  if  initially  an  atomic  O-N^  mixture 
were  used,  ^his  might  be  achieved  with  an  initial  Oj-N^  mixture.  Upon 
shock  heating,  would  quickly  dissociate  to  and  O,  thereby  preparing 

an  initial  O  concentration  before  the  N^  vibration  relaxation  has  progressed 
significantly.  Under  these  conditions,  without  coupling,  the  NO  production 
would  increase  linearly  with  time  during  N^  vibrational  relaxation.  However, 
with  coupling,  initially  the  chemical  rate  constant  would  be  less  and  the 
NO  formation  would  show  an  incubation. 

Comparison  of  Computed  and  Experimental  NO  Chemistry 

As  described  above,  the  chemical  kinetics  of  NO  production  is  a 
two-step  process  in  which  0£  is  first  dissociated,  then  the  NO  is  formed  in 
the  atom  shuttle  reactions.  The  direct  bimolecular  Reaction  (1)  is  deleted. 

In  initial  attempts  to  make  a  detailed  comparison  of  the  predicted  radiation 
with  the  experimental  data,  the  detailed  fit  was  rather  poor.  Some  prelim¬ 
inary  estimates  indicated  the  effects  of  shock  velocity  attenuation  and 
boundary  layer  on  the  shock  tube  wall  could  not  be  ignored.  A  calculating 
procedure  was  developed  which  added  the  effects  of  the  side  wall  boundary 
layer  on  the  chemistry.  Appendix  C  presents  a  detailed  derivation  of  the 
computer  program  used. 

Computer  calculations  arc  presented  in  Figs.  20,  2 1  and  22.  The 
dashed  curves  show  the  computed  times  with  the  rates  in  Table  I  (bimolccu- 
lar  reaction  omitted),  the  solid  curves  use  a  factor  of  2,5  smaller  rate 
constant  for  Reaction  (2). 
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SHOCK  VELOCITY ~  mm/fi  **c. 

Fig.  22  Time  to  3/4  peak  NO  emission  corrected  for  electronic  response 
time.  The  dashed  curve  uses  the  reactions  given  in  Table  I  with 
reaction  1  omitted.  The  solid  curve  uses  a  2.  5  times  smaller 
value  for  k,  . 
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The  major  correction  to  the  experimental  data  is  due  to  the 
finite  electronic  response  time.  The  correction  for  electronic  response 
time  becomes  significant  at  the  higher  shock  velocities  where  the  chemistry 
proceeds  rapidly.  Using  the  measured  electronic  response  time  t ■  and 
assuming  a  simple  RC  circuit  (t  =  RC),  then  the  correction  to  the  data  can 
be  obtained  from  the  relation 


Sc  =  S  +  Te(dS/dt) 


(18) 


S  is  the  observed  signal  output  on  the  oscillogram  trace,  and  Sc  is  the 
corrected  signal.  The  value  of  varied  from  0.  3  to  2  microseconds  for 
different  runs,  the  smaller  value  applying  to  the  Ce:Au  system.  The  data 
in  Figs.  20,  21  and  22  have  been  corrected  for  the  electronic  response  time. 
The  change  in  the  data  due  to  electronic  response  time  corrections  can  be 
obtained  by  comparing  Figs.  20,  21  and  22  with  Figs.  13,  14  and  15. 

In  Figs.  20,  21  and  22  the  comparison  of  the  solid  curves  with  the 
data  is  quite  good,  except  for  times  longer  than  15  microseconds.  At  these 
times  the  corrections  vo  the  clata  due  to  shock  velocity  attenuation  become 
important.  Since  we  did  not  correct  for  this  effect,  the  data  at  later  times 
were  not  used  for  analysis.  The  solid  and  the  dashed  curves  show  computed 
times  as  a  function  of  shock  velocity  and  initial  pressure  with  the  rate  con¬ 
stant  k^  as  a  parameter.  We  used  the  following  procedure  for  determin¬ 
ing  our  experimental  rate  constant  product  for  each  of  the  experimental 
points.  The  NO  concentration  with  time  is  obtained  by  integrating 


Reaction  (2) 


NO 


[  O  ]  [  N  ^  ]  dt  * 


(19) 


The  factor  of  2  is  due  to  the  additional  NO  production  through  Reaction  (3). 
N^  is  the  original  concentration,  while  O  is  produced  by  dissociation 


*  f  2t°2){k4,02[°2)  +k4,N2[N2)  +k4,o[°]}d'' 


Insetting  Eq.  (20)  into  Eq.  (19)  and  substituting  the  rate  constant  expres¬ 
sions  from  Table  I,  we  obtain 


(21) 


(NO)  =  4  [k°  k°)  f  [N2]  exp[-T2/T)  dt'  . 

.V  °  .  /T  \V2 

•  J  [02]|z[N2]  +  9[02]  +25[0])^TUr-j  exp[-T4/T]dt" 


T  is  the  temperature:  T2  =  39,  000°K  and  T4  =  59,  000°K,  corresponding 
to  the  activation  energies  for  Reactions  (2)  and  (4^  respectively. 
k°  -  7  x  1013  and  k°  =  2.  5  x  1011  cm3/mole-sec  are  the  coefficients  of 
the  rate  constant  expressions  for  Reactions  (2)  and  (4), respectively.  The 
NO  concentration  is  shown  to  be  proportional  to  the  product  [k2  k4)  in 
Eq.  (21).  For  each  shock  tube  run  we  measured  three  experimental  times 
at  the  occurrence  of  l/4,  l/2  and  3/4  peak  NO  emission.  Experimental 
values  for  [k°  k°]  were  determined  from  these  time  measurements  in  the 
following  manner.  The  solid  and  dashed  curves  in  Figs.  20.  21  °and  22 
are  the  computer  calculated  times  for  two  different  values  of  [  k2  k4] . 
These  times  were  related  by  the  expression 


(1  +  atd) 


(22) 


where  the  subscripts  s  and  d  indicate  values  for  the  solid  and  dashed 
curves  respectively.  We  then  solved  for  a  ,  the  only  unknown  in  Eq.  (22), 
as  a  function  of  shock  velocity.  Each  experimental  time  measurement  was 
substituted  in  the  expression 


(1  4  ate) 


=  [k^l  <>  +  a,d> 

a 


(23) 


where  the  subscript  e  refers  to  experimenta.  data.  For  a  measured 
experimental  time  and  shock  velocity  the  only  unknown  quantity  in  Eq.  (2  3) 

is  [k°  k°]  . 
e 

The  ratio  of  the  experimental  rate  coefficient  product  [k2.k4]  to 

e 

the  rate  coefficient  product  given  in  Table  I  [k2.k4]  is  plotted  as  a  function 

d 

of  shock  velocity  in  Fig.  2  3.  Also  shown  are  gas  temperatures,  obtained 
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Fi«z  23  The  ratio  of  the  experimental  value  of  the  rate  constant  product 
®*  r  i.o  ,.Ot  TahU  T.  The  various  gas  temperatures  are 


[kjkj]  to  that  of  Table  I. 
also  shown. 
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from  thf  computer  calculations,  at  several  NO  emission  levels.  The  solid 
line  represents  the  average  value  and  the  dashed  lines  represent  the  standard 
deviation.  A  small  number  of  the  points  shown  in  Fig.  23  were  obtained  from 
poor  data  and  were  given  a  larger  uncertainty  factor.  The  data  seems  to  have 
a  downward  trend  with  increasing  shock  velocity,  which  can  be  fit  with  a  T 
dependence.  However,  the  scatter  in  the  data  and  errors  from  fits  to  such 
a  large  exponential  temperature  dependence  are  too  large  to  be  certain  of 
this  trend.  In  addition,  this  downward  trend  could  be  due  to  shock  wave 
attenuation  effects,  which  would  increase  the  rate  constants  at  the  lower 
shjck  speeds,  where  the  event  being  measured  occurred  further  from  the 
shock  front. 

Our  experimentally  measured  rate  coefficient  product,  [k^k^]  is 

0.40  ±  0.  12  of  the  rates  given  in  Table  I.  This  product  involves  four  rate 

constant  expressions.  There  is  about  a  factor  2  uncertainty  in  the  published 

*  23  33 

rate  constant  expressions  for  k^  q  ,  k^  q  ant*  ^2  *  ^he  va*ue  ^or  ^4, 

-  2 


is  an  upper  limit. J  If  k^  ^ 


were  dropped  by  a  factor  of  3  ,  then  the  O 2 


dissociation  by  N2  would  be  negligible.  The  rate  constant  product  above 

would  decrease  by  a  factor  2,  producing  essential  agreement  between  the 

rate  shown  in  Tabie  I  and  our  experimental  results.  A  decrease  of  k2  or 

k.  «bva  factor  of  2  would  also  bring  the  rate  shown  in  Table  I  in  agree- 
4, 0  ; 

ment  with  our  experimental  measurements.  The  values  of  the  individual 
rate  constants  cannot  be  obtained  from  our  data  and  would  require  further 
measurements  with  various  N2™02  mixtures. 


Conclusion 

We  have  reported  on  infrared  measurements  of  NO  formation  in 

shock  heated  air  from  which  the  important  reaction  rate  constants  were 

evaluated.  Our  experimental  value  of  the  NO  5.  3  micron  fundamental  band 

-1  -2 

intensity  is  124  ±  22  atm  cm  . 

The  reaction  N^  +  02  —  2NO  proposed  by  Freedman  and  Daiber 
and  Wray  and  Teare^  ie  incorrect.  We  obtained  data  repeating  their 
experiment  on  shock  heated  NO  to  which  they  attributed  the  reverse  of  the 
above  reaction.  We  have  shown  that  the  NO  depletion  was  due  mainly  to 
2NO  -  N20  +  O  . 
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Without  Reaction  (1),  NO  production  is  initially  slow  until  formation 
of  atomic  O  ,  but  then  proceeds  rapidly  by  the  atom  shuttle  reactions.  The 
NO  production  depends  on  the  rate  constant  product,  k^. 

It  was  shown  that  corrections  due  to  the  shock  tube  boundary  layer 
growth  must  be  included  in  order  to  obtain  detailed  agreement  between  experi¬ 
ment  and  calculations.  The  method  is  developed  in  Appendix  C.  An  analysis 
of  the  vibrational  relaxation  of  NO-02  mixtures  is  presented  in  Appendix  A. 

NO  vibration  lag  is  small  and  dors  not  affect  the  data  analysis. 


APPENDIX  A 


NO  VIBRATION  RELAXATION  IN  AIR 

The  theory  of  vibrational  relaxation  in  a  mixture  of  diatomic  molecules 
is  reviewed.  Numerical  estimates  are  then  presented  for  the  relaxation  of 
an  N^-O^-NO  mixture.  Our  specific  problem  is  to  determine  the  vibrational 
temperature  history  of  a  trace  amount  (less  than  10"2  concentration)  of  NO 
in  temperatures  ranging  from  2000°K  to  6500°K,  corresponding  to  shock 
velocities  from  2.0  to  4.0  mm  per  microsecond,  see  Fig.  1. 

In  a  mixture  of  diatomic  gases,  the  vibrational  relaxation  rate  of 
each  component  can  differ.  This  leads  to  added  complexity  in  the  number 
of  ways  in  which  to  excite  vibration.  Not  only  is  there  the  usual  vibrational 
excitation  in  which  collisional  energy  is  converted  to  vibration,  designated 
as  a  (TV)  process,  but  vibrational  energy  can  be  transferred  from  one  gas 
species  to  another,  designated  as  a  (VV)  process.  Consider  the  example  of 
an  N2-Oz  mixture  for  which  there  are  five  vibrational  excitation  reactions. 


o2  *  o2  -  02  t  0* 

(A-l) 

N2  +  N2  -  N2  +  N* 

(A-2) 

°2  +  N2  -  °2  +  N* 

(A-3) 

°2 1  n2  -  o*  ♦  n2 

(A-4) 

°2  +  N2-°2tN2 

(A-5) 

1  \c  asterisk  indicates  vibrational  excitation.  The  first  four 

are  (TV)  pro 

cesses;  the  last  one  a  (VV)  process. 

In  a  mixture  of  the  three  molecules  NO,  N2  and  02,  there  are  12 
reaction  paths;  3  pure  gas  (TV)  processes,  similar  to  Reaction  (A-l); 

3  (VV)  processes,  like  Reaction  (A-5);  and  6  two-component  (TV)  processes 
like  Reactions  (A-3)  and  (A- 4).  Of  these,  7  have  been  experimentally  studied. 
Extensive  measurements  of  the  vibrational  relaxation  of  O  34' 35N  2^' 36and 

w  b 
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NO37,  38  have  been  reported.  Vibrational  relaxation  in  0,-N,30,  36  and 

31  2  2 
NO-N^  mixtures  has  been  observed.  The  experimental  relaxation  times 

for  these  gases  for  a  pressure  of  one  atmosphere  are  plotted  as  a  function 

of  T  3  in  Fig.  A-l  as  the  solid  lines. 

Estimates  of  the  rate  constants  for  the  other  reactions  were  obtained 
from  theory.  The  dashed  curves  in  Fig.  A-l  give  the  vibrational  relaxation 
times  at  one  atmosphere  pressure  as  a  function  of  temperature.  The 
measured  NO-N2  and  C^-N^  vibration  exchange  rates  were  1/10  the 
theoretically  predicted  rate.39  Thus,  the  02~N0  vibration  exchange  rate 
was  also  set  at  1/10  the  theoretical  rate.  The  estimated  (TV)  reaction  rates 
are  not  important  reactions  for  relaxation  in  a  trace  NO-air  mixture. 

By  combining  the  rate  constants  obtained  *Vom  experimental  measure¬ 
ments  and  theoretical  predictions,  we  can  estimate  NO  relaxation  in  air. 

We  followed  the  mathematical  treament  of  Schwartz,  Slawsky  and  Herzfeld39 
for  vibrational  relaxation  in  a  mixture  of  diatomic  molecules.  Starting  with 
Eq.  (22)  of  Ref.  39,  the  rate  of  production  of  vibrational  energy,  E  ,  in  a 
mixture  of  three  molecules  designated  by  the  subscripts  a,  b  and  c  becomes 


$ 

ip  is  the  mole  fraction  of  the  ath  molecule.  Thus,  ip  +  tp,  +  tp  =1.  E  is 
a  Ta  Tb  rc  a 

the  equilibrium  vibrational  energy  at  the  translational  temperature,  T. 

Since  a  harmonic  oscillator  potential  was  used  in  the  derivation  of  Eq.  (A-6), 


^Equation  (22)  of  Ref.  39  has  a  few  mistakes.  Their  previous  Eq.  (20)  is 
correct  and  would  give  the  formulation  above  for  their  Eq,  (22). 
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(A-7) 


E&  has  the  simple  form 


+ 

E  a 
a 


tb  Nhv 
ra  a 

Ti^/UT 

e  - 1 


N  is  the  total  number  of  molecules  per  cm  .  The  harmonic  oscillator 

molecular  potential  has  equally  spaced  vibrational  levels  with  energy  spacing 

(hv).  M..  is  the  number  of  collisions  per  second  that  an  ith  molecule  makes 

'  ij  3 

with  jth  molecules,  assuming  there  are  N  jth  molecules  per  cm  .  is 

obtained  from  the  relation 


M.j  «  N  c  tRj  +  Rj)  = 


nV 


(rn.  +  m.)  <R.  +  R.)‘ 
nij  i  J  *  J 


(A-8) 


"c  is  the  mean  particle  velocity  based  on  the  reduced  mass  of  an  i  and  j 
molecule;  m  and  R  are  the  mass  and  radius.  P.j  is  the  probability  that 
in  a  collision  of  an  ith  and  jth  molecule,  the  ith  molecule  makes  the  vibra¬ 
tional  transition  (v^  =  1  to  v.  *  0)  transferring  energy  to  the  translation  mode. 
Vj  is  the  vibrational  quantum  number.  P?j  is  the  probability  that  in  a  collision 
of  an  ith  and  jth  molecule,  vibrational  energy  is  exchanged  in  the  transition 
(v.  =  1  to  v.  =  0)  and  (v^.=  0  to  v.  =  1).  Differential  equations  for  the  relaxa¬ 
tion  of  E^  and  Ec  are  similar  to  Eq.  (A-6)  and  are  obtained  by  suitably 
changing  the  subscripts.  There  are  two  simplifications;  from  Eq.  (A-8) 

M..  =  M..  (A-9) 

ij  J» 

and  theory  predicts  ^ 

-Mv.-vJAT 

(A- 10) 


-h(v.  — v.  )/kT 
P7.  =  P7.  e  1  J 


No  general  relation  between  P.j  and  Pj.  will  be  assumed.  It  can  bo  shown 
that  a  system  with  an  harmonic  oscillator  potential  and  which  satisfies  the 
harmonic  oscillator  selection  rules  has  a  vibrational  temperature  Tv 
during  the  entire  relaxation  process,  which  is  related  to  the  vibrational 
energy 


17 


E  = 
a 


4>  Nhv 
a  a 

hv  /kT 
a  v  . 
e  -1 


(A-l  1 ) 
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In  order  to  relate  the  set  of  relaxation  equations  to  the  experimental 

measurements,  we  transfoim  equations  of  the  (A- 6)  type  to  those  in  terms 

of  the  vibrational  relaxation  times.  For  a  (TV)  process,  the  vibrational 

relaxation  time  T..  is  defined 
U 

-hy./kT 

1/tj.  =  M..P..U-C  1  )  .  (A- 12) 

£ 

For  a  (VV)  process,  the  relaxation  time  r.j  is  defined 

l/rC.  =  Ki..PC.  .  (A- 1 

Comparing  Eq.  (A- 13)  with  Eq.  (A- 10),  we  obtain 


(A- 13) 


where  9  ,  =  e 

i 


-h(i/.-y.)/kT 

t°.  =tC.  e  1  J  =  t..(G  ,/G 
ij  ji  Ji'  /  J 

-hv./kT 

.  We  introduce  the  new  dependent  variable 
X.  -  (E*-  E.)/E*  , 


(A-14) 


(A- 15) 


which  is  preferable  to  the  energy  as  it  varies  exponentially  with  time  for 
systems  with  only  (TV)  processes. 

Substituting  Eqs.  (A-12)  through  (A-15)  into  Eq.  (A-6)  and  specifying 
a  NO-N^-O^  mixture,  one  obtains  the  relaxation  equations 


dX0,  'o,  ^  NO 

- 1= - - —  + - —  + -  xQ 

dt  to2-n2  t<Vn0  2  to2-n2 


XN2(,-XOi>  X02"-XN2» 


d-e  )- 


4,  NO 

'xno<‘-xo2> 

e 

to2-no 

O 

<3> 

1 

_ 1 

02'  NO 
19  NO 


(A- 16a) 


<‘-%2> 
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dt  T 


V°2  TN2'N2  VH  T°2-N2 


^NO  ‘^c/9!^  XN2(1’5<D^)  XDz(1_XN2) 

-  +  — - /'1"0N)+ 

VN0  t02-N2  *-e02  ‘-eN2  J 


(A- 16b) 


'no  'SNc/eN2{!%<l-XNO>  W’V 


dXNO. 

^O 

2 

dt 

tno-o2 

tno-n2 

e 

TNO-N2 

1  ! 

O 

z 

0 

1 

<D 

z 

^NO 

V0o/0NO> 

XNO(l"*o) 

tno-no 

+  e 

to2-no 

>-%2 

** 

~NO 
(A  -  l6c) 


"-0NO>--5 - 

tno-n. 


n2  K(1-w  w-v 


"-eNO>- 


In  order  to  obtain  numerical  solutions,  each  T  and  9  must  be  specified. 
Numerical  values  for  T  may  be  obtained  from  Fig.  A-l.  Using  published1  ^ 
molecular  constants, 


_  -2239°K/T  = 

y  NO 


2702°K/T.  _  =  -  3353°k/T 


Wc  obtain  solutions  of  Eq.  (A- 16)  for  the  specific  case 

#0  *  °-2l>  *N  =  °'79:  *N0  =  -0001  • 

2  2 


(A- 17) 


(A- 18) 


Equation  (A-l 6)  can  be  solved  once  the  kinetic  temperature  is  given  as  a 
function  of  time.  For  the  shock  tube  experiments,  the  temperature  T  is 
obtained  by  combining  Eq.  (A- 16)  with  the  Rankine-Hugoniot  relations 
across  the  shock  for  the  conservation  of  mass,  momentum  and  energy. 
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(A- 19) 


pu  =  p1u1 

x  2  ,  2 

p  +  pu  =  pj  +  p1ul 
h  +  }  u2  =  +  |  Uj2 

For  our  experimental  conditions,  the  perfect  gas  relation  applies 

P  =  PRT  ,  (A-20) 

and  the  enthalpy  per  unit  weight  h  is  given  by 


h  =  3.  5RT  +  i  (En  +  E 
P  v  -j 


n2  +  eno)  * 


(A-2 1) 


The  vibrational  energies  are  obtained  from  solutions  of  Eq.  (A- 16), 

The  set  of  equations  (A- 16,  19,  20  and  21)  was  programmed  for  a 
7090  computer.  Figure  A-2  shows  the  behavior  of  XNr~,  and  X„  with 

time,  for  a  shock  velocity  of  2.0  mm  per  microsecond,  5  torr  initial  pres¬ 
sure,  297°K  initial  temperature,  and  0.01%  NO  in  air. 

There  arc  several  important  features  to  these  curves:  (1)  initially, 
the  02,  NO  and  N2  relax  without  (VV)  processes;  (2)  (VV)  processes  become 
important  after  sufficient  Oz  excitation;  and  (3)  finally,  the  three  gases  follow 
a  single  relaxation  time.  The  NO  quickly  relaxes  to  about  60%  of  equilibrium 
and  then  follows  the  final  relaxation  rate.  T  *0_0  being  uncertain  was 

changed  by  a  factor  of  3  and  0.  3;  this  has  only  a  small  effect  on  the  NO 
relaxation.  Figures  A- 3  and  A-4  show  the  XNQ,  XQ  and  XN  variation 

with  time  for  shoc.c  velocities  of  2.  5  and  3.  0  mm  per  microsecond  using 
the  same  initial  conditions.  These  calculations  of  the  NO  vibrational 
relaxation  time  in  shock  heated  air  are  used  to  determine  the  decrease  in 
the  infrared  emission  due  to  vibration  lag.  The  infrared  emission  is  pro¬ 
portional  to  the  vibrational  energy,  see  Eq.  (4).  Thus,  the  ratio  of  the 
emission  to  that  at  equilibrium  vibration  is 


I/loq=E/E*=  1-Xno.  (A-22) 

Values  of  I/I^  presented  in  Fig,  il  were  obtained  from  Figs.  A-2,  A-3  and 
A-4,  using  Eq.  (A-22). 
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APPENDIX  U 


CALIBRATION  OF  INFRARED  SYSTEM 

There  arc  two  parts  to  this  appendix.  First,  the  mathematics  is 
deveJoped  for  calibrating  the  infrared  system.  These  formulas  are  used 
in  Section  III  to  determine  the  NO  fundamental  band  strength.  Second,  we 
show  that  the  weak  line  approximation  is  valid  for  the  analysis  of  the 
shock  tube  data. 

The  infrared  system  calibration  is  determined  from  the  molecular 

constants  for  the  5  3  micron  NO  fundamental  vibration-rotation  band,  hi 

the  general  case,  the  radiation  depends  upon  the  gas  composition,  tempera- 

40 

ture,  optical  path  density.  See,  for  example,  the  treatments  of  Breen, 

41 

and  Carpenter  and  Franzosa.  In  this  experiment,  the  NO  concentrations  and 
partial  pressures  are  low  enough  so  that  the  effects  of  self-absorption  are 
not  important.  We  shall  show  at  the  end  of  thi6  appendix  that  the  gas  is 
optically  thin;  thus,  the  intensity  is  simply  the  sum  of  the  emission  from 
each  molecule. 

1 5 

The  notation  of  Herzberg  is  used.  The  NO  ground  electronic  state 
is  a  n3/i>  j  ^  doublet  system  with  three  branches  of  two  closely  spaced 
lines.  For  our  purposes,  this  system  can  be  approximated  with  the  formulas 
for  the  symmetric  top  with  A  -  1.  The  emission  lines  in  the  three  branches 
of  the  fundamental  band  are  produced  by  the  vibrational  transitions  v>  to  v"  = 
(v'-l)  and  the  rotation  transitions 

P  branch:  J'  to  J"  =  J'  +  1  J'=  1,2,3,...  (B-la) 

Q  b)  anch:  J'  to  J"  =  J'  J'  =  1,2,3,...  (3-lb) 

R  branch;  J'  to  J"  =  J'  -  1  J'  =  2,  3,4 _  (B-lc) 

v  and  J  are  the  quantum  numbers  for  the  vibration  and  rotation  levels, 
respectively.  The  prime  and  double  prime  superscripts  correspond  to  the 
upper  and  lower  states  in  the  transition.  The  wave  numbers  of  the  rotation 
lines  for  the  P,  Q  and  R  branches  are 
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v  *  V  ,  -  2(J'+1)B  -a  [(J'  +  l)2  -  2(J'  +  l)v'-l]  J*  >1,2.3,.,.  (B-2a) 

v  -  v  ,  -  q  ( J1  +  J*  -1)  J'  =  1,2,3,...  (B-2b) 

Q  v  •  e 

V,  =  V  ,  +  2 J’  B  -a  [ J,Z  +  2 J'v'  -  1]  J'  =  2,  3,  4,  .  .  .  (B-2c) 

R  v'  c  el 

The  subscripts  P,  Q  and  R  indicate  the  branches.  At  the  band  center  the  wave 

number  u  ,  with  upper  level  v'  is 

V  ?  2-1 
„  =w  -  2v'  «  x  +  3.25  v'  w  y  ■  |  1904.  03-27.  9  v*  -  0.  0039  v'  ]  cm  (B-3) 
v'  e  e  e  e'e 

B  and  a  are  parameters  pertaining  to  the  rotation  constants.  For  NO, 

E  =  1 . 7046  cm  'l  and  a  «  0.017P  cm"1.  Note  that  the  P  and  R  branches 
e  e 

form  a  progression  of  almost  equally  spaced  lines,  except  for  three  missing 
lines  at  the  center:  the  J'  =  0  line  in  the  P  branch  and  the  J1  =  0  and  1  lines 
in  the  R  branch.  The  O  branch  with  closely  spaced  lines  is  at  the  center 
of  the  band.  There  is  an  R  branch  band  head  for  each  value  of  v'.  For  v'  =  1 , 
the  band  head  occurs  at  2022  cm"1  for  J'  =  84.  The  P  branch  does  not  have 
a  band  head  and  expands  to  smaller  wave  numbers. 

The  intensities  of  the  rotation  lines  are 


lP  Z 


v'V*  /j,2+2, 

em  P  I 


2J'\  -hc[J'Z  +  J'-l]  [Be  -  Qe(v'4)]/kT  ,hcG  /kT 


e  v’ 


(B-4a) 


T  em 

Q  U 


<  *'-0*  *  ‘\  -hc[J'2*J'-l][Be.%(v4)]/kT  _hcC  /kT 


lv' VQ  (lV  *  l\ 


e  v' 


(B-4b) 


c.mv,,,R4  *‘\  *hc  [J'2+J'-l]  [Be-.e(v'4»)AT  ,hcC  t/kT  (B.«C) 

Ir=_°7d r  e  v 

^  9  th 

G  is  the  term  difference  (in  wave  numbers)  between  the  v  vibration  level 
v 

and  the  ground  level,  v  =  0,  ^ 

Gy  c  v[u>e-(v  +  1)  we*e  +  (v2  +  l.  5  v  +  0.75)weye]  (B-5) 

Note  G  ,  =  v,,  compare  Eq.  (B-3)  and  (B-5).  C>  and  Q  are  the  partition 

v  =  1  l  v  r 

functions  for  vibration  and  rotation.  Qf  is  a  sum  over  the  P,  Q  and  R 

branches  of  the  factors  which  contain  J'. 
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♦ 


J'=l 


yL  +  2j'  ,  iy  +  1  ,  j'^-i 

*T  -rrrrr— r»  *1 


T  TT —  7  JTJ'+T)  T  J 


-he  [J'Sj'-l]  [B(i'at>(v'  +  ]  )]/kT 


(B-6) 


_ kT 

hclBp"aP(v'  +  i  )] 


Q  is  defined 
v 


Q  * 


Z  e 


-hcG  /kT  -hey.  /kT  , 

v  ,  1  >  - 1 


~(1  -e 


(B-7) 


=  0 


40 


Breen  computed  the  matrix  elements  for  the  anharmonic  potential 

of  NO  and  finds  that  c  is  not  constant  but  ha6  a  slight  variation  with  v' . 

em 

Hi6  result,  for  the  fundamental  band,change6  cem  in  Eq.  (B-4)  to 


c*  (y  ,  Jv  ,)*■--  c*  /[1-0.  01  5(v’ -1 )]  2 
l  v'  =1 '  v'  em'  1  ,J 


etn 


(B-8) 


where  c  is  the  value  of  c  for  v'  =1.  Thi6  corrected  matrix  element 
em  em 

was  not  used  in  our  calculation  as  it  has  a  negligible  effect.  At  a  tempera¬ 
ture  of  6000°K,  it  would  increase  the  total  emission  by  about  9%.  At  lower 
temperatures,  the  increase  is  even  smaller. 

Except  for  c  ,  numerica.  values  of  all  the  quantities  in  Eq.  (B-4) 
em 

are  now  specified.  Our  experimental  data  will  be  used  to  determine  c  _ . 

r  em 

We  will  use  the  following  procedure. 

In  this  experiment,  we  U6ed  broad  band  infrared  (InSb  and  Ge:Au) 
detectors  with  interference  filters  to  isolate  the  NO  band,  (see  Fig.  5). 

Let  r(y)  be  the  transmission  of  the  filter,  and  o(y)  the  response  of  the 
detector;  both  quantities  are  functions  of  wave  number  y  .  The  signal  S 
observed  by  the  detector  i6  determined  by  the  emission  of  the  gas  folded 
into  the  .response  of  the  detector  aftd  the  filter  transmission; 


S  =  v  [  or  Ip 

v'=lJ'=I 


+  otIq  +  cttIr] 


fl  is  the  solid  angle  of  the  gas  as  seen  by  the  detector.  Inserting  the 

values  of  I_,  I0  and  I_  from  Eq.  (B-4)  into  (B-9),  one  obtains 
r  U  R 


(B-9) 
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is  the  solid  angle  of  the  optical  system  using  the  full  tt-incn  cnamcii 
ror.  X  is  the  fraction  of  the  mirror  uncovered;  that  is,  V  =  «  0X0  1S 
solid  angle  actually  used.  N  is  the  NO  number  density,  and  i  3.81  c 
he  diameter  of  the  shock  tube.  The  temperature  dependence  of  the  calibi 
i  shown  in  Fig.  7  was  computed  from  the  bracketed  portion  of  Eq.  (B-10) 
height  of  the  curve  was  adjusted  for  the  best  fit  to  the  data.  The  depend- 

e  of  a  and  t  on  wave  number  is  presented  in  Fig.  5. 

TV.®  in fm r flvstem  was  also  calibrated  with  a  black  body  source  at 


I 


Validity  of  Weak  Line  Approximation 

We  will  prove  that  the  gas  is  optically  thin  over  our  experimental 
conditions  by  showing  that  the  gas  emissivity  at  the  center  of  each  spectral 
line  is  much  less  than  unity.  The  emissivity,  £  ,  of  a  layer  of  radiating 
gas  can  be  written 

(  :  1  •  exp  (-a  Nf  p)(l  -  exp  [  -hcv/kT] )  (B-13) 

p  is  the  absorption  coefficient  per  particle  and  Nf  is  the  thickness  of  the 
gas  sample  in  particles  per  cm2;  a  is  a  geometric  coefficient  of  the  order 
of  1.  The  factor  (1  -  exp  hcy/kT)  varies  between  .  4  and  .  9  for  our  experi¬ 
mental  conditions.  For  €  «  1.  Eq.  (B-13)  reduces  to 

€  «  N1  n  (B-14) 

Within  the  accuracy  of  this  derivation,  the  factor  a(l -exp[ he v/WT] )  is  the 
order  of  unity  and  has  been  dropped.  The  emissivity  for  the  absorbing 
state  (v",  J")  is 


«  =  N*  M(v„,  Jtl)=  H  *P(V  +  l)(273°K/T)CabBbo^R 


■  (Gvll+[  Be-  Qe(v"+| )]  [  J"  +J"-1  ]  )hc  v/kT 


Q  Q 

r  v 


(B-l  5) 


The  factor  j  is  due  to  the  doublet  system.  An  R  branch  line  is  used  as 
it  has  the  maximum  absorption.  P  is  the  pressure;  P  is  the  partial 
pressure  of  NO;  bQ  is  the  line  width  factor.  Cabg  is  determined  from  the 
measured  integral  absorption  over  the  band. 


w  w 

EE  / 

v"  0  J=1  band 


V  J"d>'«cabs,'l  =  124  atm^cm"2  (B-16) 


P..Q.R 

-1 


Using  =  1876  cm”  ,  the  band  center  wave  number,  Cabg  becomes 


124 


=  0.  066  atm  'em  * 


abs  -  rrm) 

Let  b{v)  be  the  spectral  shape  of  the  line.  In  the  general  case, 
there  is  both  collision  and  Doppler  broadening,  the  line  shape  becomes 


(B- 17) 
where 


-64- 


u/  \  —  ^  f  Yl^,  /  Q.  ^ 

b(v)-  7 t/t 


exp  [  -  (4  £n2)(  v  '  -  vR)2/p2  ) 
(q/2)2  +  (v-v  ')2 


dv'  (B- 18) 


a  and  p  are  the  full  widths  of  the  line  at  half  amplitude  for  collision  and 

Doppler  broadening  respectively.  Equation  (B-18)  is  normalized  for  each 

line  so  that  b(v)  dv  =  1.  At  the  center  of  the  line,  say  b  ,  where  the 

42  ° 

absorption  is  maximum,  Eq.  (B-18)  becomes 

bQ  =£  b(v=vR)  =  <~)  exp [  (-^-)  (£n2)]  erfc  (B-19) 

The  collision  broadened  line  width  at  standard  temperature  and  density,  a  , 

- 1  ° 

was  measured  for  collisions  with  Ar;  a  =  0.  06  cm  .  We  assume  the  same 

o 

value  of  qq  for  collisions  with  and  The  dependence  of  a  on 

temperature  and  total  pressure  is 


P 

a  =  p—  a 
o 


,ilF 


(B-20) 


where  Tq  =  273  K  and  Pq  =  1  atmosphere  pressure.  The  Doppler  half  width 


P  =  (.0040  cm"1)  (vr/vj)  ^T/Tq 


(B-21) 


For  our  experimental  conditions,  the  line  width  is  less  than  0.  025  cm 
This  should  be  compared  to  the  average  spacing  between  lines  of  3.  4  cm 
Thus,  the  spacing  is  100  times  the  line  widths;  the  overlap  of  lines  is  negli¬ 
gible.  Also  negligible  is  the  overlapping  of  lines  of  the  P  and  R  progressions 
for  different  vibration  levels,  v". 

By  optimizing  Eq.  (B-15),  one  can  show  that  the  e  has  the  maximum 
value  for  the  R  branch  line  with  v"  =  0  and  J"  =  Qr/2.  Using  these  values 

of  J"  and  v",  and  combining  Eqs.  (B-16)  and  (B-15),  one  obtains  the  desired 
form  for  the  maximum  emissivity. 


€  =  62  £  b 


fTo  hc  Be 

iT  “FT" 


r* 


-hcv, /kT 
(1-e  ) 


(B-22) 


- '  Ka 


When  evaluating  Eqs.  (B-19)  and  (B-22)  over  the  experimental  conditions 
shown  in  Pigs.  1,  2  and  3,  it  is  found  that  e  <  0.  1.  Thus,  the  gas  is 
optically  thin  at  the  center  of  each  line. 


APPENDIX  C 


EFFECTS  OF  BOUNDARY  LAYER  GROWTH  ON  SHOCK  TUBE  WALLS 

Many  investigations  of  boundary  layer  growth  on  shock  tube  walls 

have  been  reported.  These  have  mainly  been  concerned  with  the  effect  of 

the  boundary  layer  on  shock  tube  test  times  and  on  the  variation  of  equilibrium 

conditions  along  the  shock  tube  axis  caused  by  the  boundary  layer  growth. 

43 

Mirels  discussed  these  effects  in  terms  of  a  "limiting  separation  distance,  " 

44 

first  observed  by  Duff.  This  concept  is  depicted  in  Fig.  C-l,  taken  from 
Ref.  43.  A  normal  shock  wave  with  velocity  Ug  is  advancing  into  stationary 
gas;  the  profiles  in  the  figure  show  velocities  relative  to  the  shock  front, 
and  in  this  coordinate  system  the  shock  tube  walls  are  moving  with  the  shock 
velocity.  Ahead  of  the  shock  the  relative  velocity  of  the  gas  is  supersonic, 
while  immediately  behind  the  shock  the  velocity,  ucq  ,  is  given  by  Peoueo  = 
PjUg.  As  the  boundary  layer  develops,  u,  and  pe  change,  such  that  the 
integral  of  pu  over  the  entire  shock  tube  area  is  conserved.  However,  in 
this  coordinate  system  the  boundary  layer  is  a  region  in  which  p  >  pe  and 
u  >  uc  ,  corresponding  to  a  negative  displacement  thickness.  Thus  the  gas 
which  has  not  yet  entered  the  boundary  layer  behaves  as  if  it  were  in  a 
divergent  channel,  resulting  in  a  decrease  of  ug  with  increasing  distance  * 
behind  the  shock,  since  ueQ  is  necessarily  subsonic.  Ultimately  the  bound¬ 
ary  layer  grows  such  that  the  integral  of  pu  across  the  boundary  layer 
area  is  equal  to  the  integral  of  Peoueo  across  the  entire  tube  area.  At 
this  point,  a  distance  x^  behind  the  shock  wave,  the  central  test  gas  is 
stagnated. 

In  terms  of  the  operation  of  a  shock  tube,  the  phenomena  may  be 
described  as  follows:  after  rupture  of  the  diaphragm  and  establishment  of 
a  plane  shock  wave,  the  distance  between  the  shock  and  the  contact  surface 
initially  grows  with  time.  Once  this  distance  has  reached  the  limiting 
separation  distance  x £  ,  the  rate  of  flow  of  cold  gas  through  the  fchock 
wave  is  exactly  balanced  by  the  rate  of  outflow  through  the  boundary  layer 
at  x  =  X£  ,  and  no  further  increase  of  length  of  the  "test  slug"  occurs. 

When  a  shock  tube  of  conventional  length  is  operated  at  sufficiently  high 


-67- 


BOUNDARY 

LAYER 


CONTACT 

SURFACE 


SHOCK 


Fig.  C-l  A  schematic  representation  of  the  quantities  used  in  calculating 
the  effects  of  the  side  wail  boundary  layer. 
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pressure  this  maximum  length  of  test  slug  is  not  achieved,  and  the 

variation  of  u  with  distance  behind  the  shock  may  be  negligible  during  the 
e 

experimental  test  time.  However,  for  many  chemical  kinetics  experiments 

the  shock  tube  initial  pressure  must  be  kept  low  in  order  to  obtain  a  finite 

relaxation  zone.  Under  such  conditions  the  experimental  observations  may 

continue  over  almost  the  entire  available  test  time,  and  significant  errors 

may  then  be  incurred  in  the  data  analysis  if  the  variation  of  ue  is  neglected. 

The  variation  of  u  has  a  two-fold  effect  on  the  interpretation  of 
e 

shock  tube  kinetics  data.  In  the  first  place  the  deceleration  of  the  subsonic 
flow  relative  to  the  shock  wave  increases  the  static  enthalpy  of  the  test  gas, 
and  this  can  have  a  marked  effect  on  the  reaction  rates  of  reactions  with 
steep  temperature  dependence.  Secondly,  the  deceleration  results  in 
increased  flow  time  for  a  given  gas  sample  as  it  flows  from  the  shock  wave 
to  the  point  of  observation.  The  distance  from  the  shock  to  the  point  of 
observation  is 


x 


9 


where  tp  is  the  particle  time  measured  from  the  shock  front.  In  the 
laboratory  coordinate  system,  the  time  tj^  is  measured  from  the  arrival 
of  the  shock  at  the  observation  station,  so  that 


In  the  absence  of  boundary  layer  effects,  the  ratio  ug/ Ug  =  Pj/Pe»  ar*d  is 

affected  only  by  changes  in  p  as  chemical  reactions  proceed,  but  when 

e 

ug  can  decrease  towards  zero  it  is  clear  that  dtp/dt^  can  become  very 
large  as  x  approaches  the  value  x^  . 

In  the  present  study  of  NO  formation  in  shock  heated  air  it  was 
apparent  that  some  of  the  data  obtained  at  low  initial  pressures  were 
affected  by  the  boundary  layer  growth.  In  comparing  the  data  with 
theoretical  calculations,  it  was  thus  necessary  to  include  the  effective 
area  variation  of  the  channel  in  which  the  reacting  gas  flows  after  passage 
through  the  shock  wave.  In  the  conventional  constant  area  1 -dimensional 
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normal  shock  chemistry  calculation  it  is  usual"  to  integrate  the  various 
chemical  rate  equations  subject  to  constraints  imposed  by  the  shock 
conservation  equations,  which  can  be  written  in  the  form 

Pe  =  Pi  +P  1  ^'Pl/Pe* 


and 


h  =  h,  +  1  U  ‘ 
e  *  2  ■ 


I1  - 


since  pe  ue  =  Pl  Ug 


In  the  present  case,  it  was  more  expedient  to  modify  a  stream-tube 
program  to  permit  integration  of  the  chemical  rate  equations  subject  to 
constraints  imposed  by  a  specified  variation  of  area  with  distance. 

Thus  the  governing  equations  become 

h  ♦  1 

e 


ana 


with 


dPe 

-HtT 


:  *  hi  +  \  v< 2 

(C-l) 

-Pi  °.Ao 

(C-2) 

due 

“Peue*3t^  ’ 

(C-3) 

p  RT 

(C-4) 

Equations  C-l,  C-2  and  C-4  car.  be  written  in  the  form 


dh 


and 


=  -u 

dPe 

+  j_ 

du 

e 

u 

e 

dPc 

1 

dze 

“3t^ 

ze 

dtP 

due 
e  Ht^ 


,  1  dA.  _  o 

.  1  1 

P~  +  T~ 
He  P  e 


dT 

< 

■an: 


(C-5) 


(C-6) 


(0-7) 


-70- 


Eliminating  the  derivatives  of  ue*  pg  and  pe  from  Eqs.  C-3,  C- 5,  0.-6 
and  C-7,  we  obtain 


dh  (  . 
ell 

Ttr  “T 

P  (  u 


i  dz  , 

1  e  ,  1 

+  -  -T7 -  +  TS - 

ze  ctP  Te 


dT 


7T  0  (c'8) 


Now  for  any  given  set  of  reactions  and  reaction  rates  the  enthalpy  derivative 

dh  dl 

can  be  expressed  in  the  form  =  BL  +  B.,  ,  where  the  c°efficlents 

B  and  B,  involve  derivatives  of  species  concentrations  -*r-i  and  of  vibra- 

1  l  dT  P 

tional  temperature  -rj — ,  and  these  are  all  functions  of  local  gas  proper- 

dtP  dz  \dX. 

ties  only.  The  value  of  ^  is  also  3  function  of  local  gaS  pr0p" 


ertie6  only. 

Consequency,  Eq.  (8)  provides  a  direct  specification 
provided  —  is  known.  Values  of  the  derivatives  of  he. 


and 


pe  then  follow  directly  from  the  above  equations. 

The  Avco  Everett  Research  Laboratory  stream-tube  program  was 
modified  along  the  lines  indicated  above,  with  the  area  variation  of  the 
diverging  channel  specified  with  the  aid  of  Mirels'  work.  For  a  laminar 
boundary  layer  and  a  given  value  of  limiting  separation  distance  x^  , 
the  area  variation  experienced  by  the  shocked  gas  is  given  by 


A  /A  =  (p  u  )/(p  U  )  =  1  -  (x/x  -  )n,  with  n  =  0.  5.  (C-9) 

The  same  expression  with  n  =  0.  8  is  applicable  for  a  turbulent  boundary 
layer  provided  an  appropriate  value  of  X£  is  used. 

Only  one  difficulty  remains  in  the  application  of  this  model  for 
interpretation  of  the  present  experiments.  It  may  be  noted  that 

1  dA  ue  dA  A 

=  -TT-XT  -TT 

tends  to  a  value  of  irfinity  when  x  approaches  zero.  The  resulting 
computational  problems  can  be  circumvented  by  assuming  that  in  the 
region  0  £  x  <  xm  the  area  variation  is  linear,  with  A/Ao  =  1  +  Kx,  where 
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K  is  determined  from|l  +  Kxmj  -l  -  (xm^x£^  ‘  *n  choo6'n8  an 

appropriate  value  of  xm,  some  recourse  may  De  made  to  the  physical 
situation  in  the  experiment.  The  theoretical  analysis  of  the  boundary 
layer  growth  assumes  that  this  growth  commences  immediately  behind 
a  plane  shock  wave  of  infinitesimal  thickness,  whereas  the  experiment 
involves  a  slightly  curved  shock  wave  for  which-the  translatioival/rotational 
relaxation  thickness  is  approximately  four  times  the  mean  free  path  in 
the  unshocked  gas. 

For  the  1.  5  inch  diameter  shock  tube  the  translational/rotational 
relaxation  thickness  is  6g  »  0.02/pj  cm,  with  pl  *.n  torr.  For  the 
curvature,  the  effective  thickness  is  •  0.  1  / yfT[  cm,  again  with  Pj  in 
torr.  It  seems  physically  reasonable  to  choose  xm  of  the  same  order  of 
magnitude  as  either  of  these  experimental  quantities,  and  the  present 
calculations  use  a  geometric  mean, 

xm  =^5c  6s  *  '  045/Pl  °*  75  cm* 

Check  calculations  were  made  to  ensure  that  the  results  were  in  no  way 
sensitive  to  this  choice. 

Several  kinetics  calculations  were  made  for  comparison  with  the  data 
of  the  present  experiments,  using  Eqs.  (C-8)  and  (C-9).  The  values  of  x^ 
used  are  shown  as  functions  of  shock  speed  in  Fig.  C-2.  The  solid  curves 
A  and  B  are  obtained  from  Mirels’  theory43  for  a  shock  tube  of  1 . 5  inches 
diameter  with  a  laminar  boundary  layer,  for  initial  pressures  of  1  torr  and 
5  torr  respectively.  a  given  shock  speed  the  theoretical  value  of  X£ 
scales  linearly  with  Pl  .  However,  at  an  observation  station  18  feet  from 
the  shock  tube  diaphragm,  (as  in  the  present  experiments)  the  theoretical 
length  of  the  "test  slug"  would  be  approximately  40%  of  X£  for  Pj  =  5  torr, 
and  almost  100%  of  x*  for  P[  =  1  torr.  The  prediction  for  the  1  torr  case 
is  in  quite  good  agreement  with  the  experimental  test  times,  which  are  shown 
(converted  into  test  lengths)  as  the  curve  C  of  Fig.  C-l.  The  effect  of 
using  curve  A  or  curve  C  of  Fig.  C-2  in  the  calculations  is  shown  in 
Fig.  C-3.  For  pj  =  5  torr,  the  predictions  using  a  laminar  and  a  turbulent 
boundary  layer  were  not  in  good  agreement  with  the  experimental  results. 
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LIMITING  SEPARATION  DISTANCE  (X«)~cm 


SHOCK  VELOCITY  ~  mm />*  sec 


Fig.  C-2  The  limiting  separation  distance,  experimental  and  theoretically 
as  a  function  of  shock. 
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C-3  A  comparison  of  calculating  the  times  to  peak  using  various 
limiting  separation  distances  as  a  function  of  shock  velocity 
The  times  to  peak,  using  no  boundary  layer  coirection,  is 
also  shown. 
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Side  wall  heat  transfer  gauge  measurements  indicated  that  the  boundary 
layer  might  be  in  a  transition  mode.  The  experimental  test  time  for 
5  tori*  ib  shown  as  curve  D  on  Fig.  C-2.  One  should  note  that  at  the  lower 
shock  velocities  the  test  time  agrees  more  closely  with  the  turbulent 
boundary  layer  predictions,  "  Fig.  C-2,  curve  B1,  but  at  higher  shock 

43 

velocities  the  data  are  in  better  agreement  with  the  laminar  prediction, 
cui  ve  B  .  The  choice  of  an  appropriate  value  of  x^  for  p^  =  5  torr  is 
very  difficult.  Calculations  were  made  using  the  laminar  and  turbulent  pre¬ 
dictions.  The  true  value  of  x^  lies  between  these  extremes,  with  the 
initial  variation  of  area  corresponding  to  a  laminar  boundary  layer  growth 
followed  by  a  transition  to  a  turbulent  boundary  layer  at  some  value  of 
x  <  Xg  .  The  effect  of  using  either  the  laminar  or  turbulent  boundary  layer 
on  the  calculations  is  shown  in  Fig.  C-3  by  the  curves  labeled  B  and  B*. 

All  the  calculations  discussed  below  use  tha  rate  constants  shown  in 

Table  I,  except  for  modification  of  the  catalytic  efficiencies  of  0?  and  O 

47  C 

in  dissociating  O ^  and  deletion  of  the  bimolecular  reaction  +  O ^  — •  2NO. 

Figure  C-3  shows  the  results  of  several  calculations  of  time  to 
peak  infrared  emission  from  the  NO  formed  behind  the  normal  shock  wave 
in  the  presen*  experiments.  The  calculations  are  for  =  1  torr  and 
5  torr.  The  dotted-^iFhed  curves  are  for  one-dimensional  constant  area 
flow,  with  no  boundary  layer  correction.  These  two  curves  are  separated 
by  a  factor  of  5  over  the  entire  range  plotted,  since  binary  reactions  pre¬ 
dominate  (at  the  low  speed  end  of  the  pj  =  5  torr  curve  the  3  body  reactions 
are  just  becoming  significant  near  the  time  of  peak  infrared  emission)  in 
determining  the  gas  behavior,  resulting  in  flow  times  which  scale  inversely 
with  Pj  for  a  given  shock  speed.  However,  when  area  variation  is  intro¬ 
duced  into  Eq.  C-8,  it  is  clear  that  this  binary  scaling  will  be  invalidated 

unless  ^  scales  inversely  with  pj ;  i.  e.  ,  unless  x*  «  p^.  For  the 
P  * 

calculations  with  area  variation  due  to  boundary  layer  growth  shown  in 

Figs.  20-22  (solid  and  dashed  curves),  the  values  of  used  do  not  satisfy 
this  scaling  relationship.  The  solid  curves  obtained  show  a  trend  which  is 
consistent  with  the  experimental  data,  and  give  reasonable  agreement  with 
those  data  for  times  to  a  fraction  of  the  peak  radiation  less  than  15  micro¬ 
seconds.  Unfortunately,  most  of  the  data  points  yielding  times  in  excess  of 
15  microseconds  are  influenced  to  Borne  extent  by  shock  speed  attenuation, 
so  that  the  boundary  layer  growth  effects  alone  cannot  provide  adequate  c  r- 
rection  over  the  entire  shock  velocity  range. 
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The  rate  of  production  of  nitric  oxide  in  shock -heated  air  was  measured  over  the  temperature  range  from 
2300° K  to  4500  °K.  The  NO  concentration  was  determined  from  the  NO  infrared  emission  in  the  5.3  micron 
fundamental  vibration- rotation  bend.  The  main  production  of  nitric  oxide  was  found  to  be  by  the  dissociation 
reaction  O2+M^420  +  M,  followed  by  the  atomic  shuttle  reactions  0+N_c2  NO  +N  andN+O/-^  NO+O. 
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The  overall  rate  constant  k  =  k2  A4  for  NO  production  was  found  to  be  0.40 ±  0.12  of  the  rate  constant  com¬ 
piled  by  Tray. 

The  rate  constant  for  the  bimolecular  reaction  N2+  Of*  2N0  was  found  to  be  less  than  one -tenth  the 
value  suggested  by  Freedman  and  Daiber.  In  order  to  explain  the  disagreement  in  the  rate  constant,  experi¬ 
ments  were  performed  in  shock-heated  nitric  oxide  where  the  depletion  of  nitric  oxide  and  the  reaction 
products  were  observed.  It  was  found  that  the  NO  depletion  occurs  mainly  via  the  reaction  2NO->N2  O+O. 

An  analysis  of  the  vibrational  relaxation  of  N0-02-N2  mixtures  is  presented  and  it  is  shown  that  NO  vibra¬ 
tion  lag  is  small  and  does  not  affect  the  data  analysis.  We  determined  that  the  integrated  band  intensity  of  the 
NO  fundamental  is  124  ±  22  atm*1  cm*2  at  S.  T.  P. 

It  is  also  shown  that  corrections  due  to  the  shock  tube  boundary  layer  must  be  included  in  order  to  obtain 
detailed  agreement  between  experiment  and  calculations.  The  procedure  for  correcting  chemical  measurements 
for  shock  tube  boundary  layer  growth  is  outlined. 
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The  rate  of  production  of  nitric  oxide  in  shock -heated  air  was  measured  over  the  temperature  range  from 
2300°K  to  4500° K.  The  NO  concentration  was  determined  from  the  NO  infrared  emission  in  the  5.3  micron 
fundamental  vibration  -  rotation  band.  The  main  production  of  nitric  oxide.was  found  to  be  by  the  dissociation 
reaction  O2+  M^*20+M,  fr ilowed  by  the  atomic  shuttle  reactions  0  +  nJ>?  NO+N  and  N+02~>N0+0. 

The  overall  rate  constant  k  =  k2  •  k^  for  NO  production  was  found  to  be  0.40  £  0.12  of  the  rate  constant  com¬ 
piled  by  Wray, 

The  rate  constant  for  the  bimolecular  reaction  N2+  O2+  2NO  was  found  to  be  less  than  one-tenth  the 
value  suggested  by  Freedman  and  Daiber,  In  order  to  explain  the  disagreement  in  the  rate  constant,  experi¬ 
ments  were  performed  in  shock -heated  nitric  oxide  where  the  depletion  of  nitric  oxide  and  the  reaction 
products  were  observed.  It  was  found  that  the  NO  depletion  occurs  mainly  via  the  reaction  2NO->  N2  04-  0. 

An  analysis  of  the  vibrational  relaxation  of  NO-O2  -N2  mixtures  is  presented  and  it  is  shorn  that  NO  vibra¬ 
tion  lag  is  small  and  does  not  affect  the  data  analysis.  We  determined  that  the  integrated  band  intensity  of  the 
NO  fundamental  is  124±  22  atm*  I  cm *2  at  S.  T.  P. 

It  is  also  shown  that  corrections  due  to  the  shock  tube  boundary  layer  must  be  included  in  order  to  obtain 
detailed  agreement  between  experiment  and  calculations.  The  procedure  for  correcting  chemical  measurements 
for  shock  tube  boundary  layer  growth  is  outlined. 
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